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Some aspects of the synchronization in coupled maps 
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Through numerical simulations we analyze the synchronization time and the Lyapunov dimension 
of a coupled map lattice consisting of a chain of chaotic logistic maps exhibiting power law interac- 
tions. From the observed behaviors we find a lower bound for the size N of the lattice, independent 
of the range and strength of the interaction, which imposes a practical lower bound in numerical sim- 
ulations for the system to be considered in the thermodynamic limit. We also observe the existence 
of a strong correlation between the averaged synchronization time and the Lyapunov dimension. 
This is an interesting result because it allows an analytical estimation of the synchronization time, 
which otherwise requires numerical simulations. 
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Coupled Map Lattices (CML's) were introduced in the 
literature as suitable models to study spatiotemporal be- 
havior of spatially extended dynamical systems. Funda- 
mentally they are systems defined on a discrete space- 
time and possessing continuous state variables. In the 
last two decades such models have received a great and 
increasing deal of interest, being applied to several non- 
linear phenomena including systems as diverse as phys- 
ical, chemical and biological lj. CML's share with real 
complex systems one of their most intriguing behavior, 
which is the possibility of synchronization. Such a phe- 
nomenon can be observed in a great variety of real sys- 
tems, going from electronic circuits until physiological 
processes, for example 0,0. 

The model we are concerned here consists of a chain 
of N identical chaotic logistic maps |^|, each one located 
at a definite site in a discrete space, and coupled be- 
tween themselves through a power law interaction 
Our aim is to go a step further in the understanding the 
synchronization of this system through the analysis of 
the synchronization time and the Lyapunov dimension. 

The model is defined as follows. At the (discrete) in- 
stant of time n the state, or amplitude, of the map lo- 
cated at the site i (i = 1,2, ■•■ ,N) is denoted by the 
continuous variable x„ . The state of the whole lat- 
tice at time n will be given by the iV-dimensional vector 
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The time evolution of the 



system is given, in matrix form, by the following map- 
ping 

x„+i = F(x„)+IF(x„), 

where F(x) = (/(x«), /(x( 2 )), . . . , /(xW)) , and I is an 
N x N matrix which specifies the coupling among the 
elemental maps. The function /(x) characterizes the ele- 
mental map, which we assume here to be the fully chaotic 
logistic one /(x) — 4x(l — x), with x € [0, 1]. We always 



assume odd N and periodic boundary conditions. The 
initial state of each elemental map is randomly chosen. 
The power law coupling is given by 
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- (1 - Sij) - r]5 i: 



where Iij arc the matrix elements of I and nj 
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the "distance" between sites i and j. The parameters e 
and a give, respectively, the strength and the range of the 
interaction and r\ = 2 £^ r~ a is a normalization factor, 
with N' = (N — l)/2. The parameter a can assume any 
value in the interval [0, oo). The extreme cases a = and 
a — > oo correspond, respectively, to global (mean field) 
and local (first neighbors) couplings. Here we restrict the 
parameter e to the interval [0, 1] in order to get all the 
individual map amplitudes into the interval [0,1]. In a 
recent work, including one of us, it was proposed that 
the parameter e could assume any nonnegative value [|| . 
Although both systems exhibit the same behavior when 
< e < 1, our choice implies that the nonlinearity of the 
model arises solely from the nonlinearity of the elemen- 
tal maps, while the choice in ref. |fj allows, for e > 1, 
also nonlinearities arising from the mod 1 operation in the 
coupling scheme. 

Among the various characterizations of synchroniza- 
tion 0] we choose the following. A system is said 
to be in a completely synchronized state at time n if 
all the elemental maps have the same amplitude, i.e., 
x„ = x^ = ■ ■ ■ = xi^ — x* . The subspace S of all these 
states (the diagonal of the whole state space of the sys- 
tem) will be called invariant if x m S S implies that x n G S 
for all n> m. In this case a synchronized regime starts 
when the system state is put on S, which is then called 
the synchronization subspace. Given specific initial con- 
ditions, the minimum value of m for which the system 
goes into this subspace is identified as the synchroniza- 
tion time t s . The necessary and sufficient condition for 
S to be invariant is that the sum Y]^ iy be independent 
of i, i.e., the sum of all elements of each line of the ma- 
trix I must be the same for all lines 0- Besides that, in 
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order to hold x^ +1 in the interval [0, 1], for all n and i, 
it is necessary that this sum be restricted to the interval 
[—1,0]. In our case Ylj lij = f° r au * an d S turns out 
to be an invariant subspace. 

If initially the lattice is not synchronized, it only will 
synchronize after some time m if the parameters a and 
e assume values into a suitable domain, for each TV, in 
the parameter space |g, |9( . The boundaries of this do- 
main can be calculated analytically from the condition 
A2 = 0, where A2 is the second largest Lyapunov ex- 
ponent of the system, which is the largest Lyapunov 
exponent transversal to the synchronization subspace 
|10| . These boundaries can be decomposed into an up- 
per line, given by e c (a, TV) = 3/2(1 — M^/ry) -1 , and 
into a lower one, given by e c (a,N) = 1/2(1 — b^/rj) -1 , 

where = 2 YZ=i cos(2irkm/N) m - a (1 < fc < TV) are 
the eigenvalues of the circulant matrix B = (rje^I—rjel) 
0,0. For instance, the domains for TV = 5 and TV = 385 
are shown in Figure ^ The shown lines are only the 
lower ones, for both iV; the upper lines do not appear 
because in the range of figure all values of e c (a,N) are 
greater than 1.0, which are outside the range of this pa- 
rameter. We can observe two distinct regions for each TV. 
For TV — 385, and if it does not start synchronized, the 
lattice synchronizes after some time only in Region (A), 
while for TV = 5 it synchronizes only in Regions (A) and 
(B). 
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FIG. 2: Averaged synchronization time <t s > versus e, with 
a = 0.0 and a = 0.6 . (a) N=5; (b) N=385. 



(t s ) versus the strength parameter s, for two values of TV 
and a 14]. For TV = 5 (Figure H^a)), and for both the 
values of a, we observe an atypical behavior for the syn- 
chronization time when e > 0.8. It would be expected that 
(t s ) were always a decreasing function of the interaction 
strength, because greater values of this parameter would 
tend to facilitate synchronization. We have constructed 
similar plots for various values of TV and a and our re- 
sults indicate that this behavior is characteristic for small 
lattices, being therefore a finite size effect of the system. 
Such a behavior was not observed with respect to the a 
parameter. From our simulations we observed that the 
"turning point" e = e^ x \ above which (t s ) starts to in- 
crease, depends only on TV (of course, with a inside the 
synchronization domain) and it can be clearly identified 
as 



4(0,TV)-e c (0,TV) 



(1) 



It turns out that the values of TV for which this atypical 
behavior disappears must correspond to > 1. From 
the above expression (or from the numerical simulations) 
we can observe that this will be satisfied for N> Nmin = 



385. The behavior of (t s ) for TV = 
12b). 



: N m i n is shown in Figure 



FIG. 1: Domains of synchronization in the parameters space. 
For TV = 5 the CML synchronizes in Regions (A) and (B) and 
does not in Region (C). For TV = 385 it synchronizes in Region 
(A) and does not in Regions (B) and (C). 

On the state space of the system it can be defined the 
quantity d n — a n y/~N, where a n is the standard devia- 
tion of the x„ around the mean (x$ ) , with the averages 
taken with respect to the index i. It turns out that this 
quantity corresponds to the distance, in the state space, 
from the point x ra to the synchronization subspace S |ll| . 
We thus use the condition d n = as a diagnostic for the 
synchronization regime. Accordingly, the synchroniza- 
tion time t s is identified with the lowest n for which d n = 

m 

Figure[3is a plot of the averaged synchronization time 
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FIG. 3: (a) (t s ) versus TV, for a = 0.0 and e = 1.0. (b) 
Lyapunov dimension D versus TV, for the same values of a 
and e. The vertical line corresponds to TV = 385. 

In Figure |3| (a) we plot the averaged synchronization 
time with varying TV with both e and a fixed. The results 
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indicate that (t s ) tends to saturate for large N. More- 
over, the saturated time does not differs significantly 
from that corresponding to TV = 385. 

Now we consider the CML in an outer vicinity of the 
synchronization domain in parameter space. The system 
does not more attain the synchronized regime if it not 
started synchronized. Figure 0Ja) is a typical plot for a 
time series of the distances d n in such a case. To make 
easier the visualization we presented our results in terms 
of y n = — logio d n - Figure 0Jb) shows the corresponding 
statistical distributions of y n . Again the results suggest 
that, for large lattices, such a distribution is practically 
independent of N and it is very well approximated by 
the distribution corresponding to iV = 385. 

All the results presented so far indicate that iV = 385, 
apart from eliminating the just mentioned finite size qual- 
itative effects, also lead to a reasonable quantitative ap- 
proximation to the behavior of the system for larger N 
values, at least in what concern the synchronization time 
behavior and the statistical distribution of the distances 
d n . So, at least in what concerns these two aspects, we 
claim that N m i n sets a practical lower bound in numerical 
simulations for the system to be considered at the ther- 
modynamic limit. With this claim we mean that both 
the qualitative and quantitative behaviors of the system 
at the thermodynamic limit N — > oo can be reasonably 
well approximated by the corresponding behavior of the 
system with N = N min . 
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FIG. 4: (a) Time series for y n with iV = 385, a = 0.6 + 10~ J 
and e = e c (0.6, 1001). (b) Distributions of y for N — 5 (open 
circles), iV = 385 (filled circles), iV = 501 (open squares), and 
N = 1001 (filled triangles), with a = 0.6 + 10~ 5 and e = 
e e (0.6,iV). 

Now, we recall that the time oscillation of d n in Figure 
EJa) is due to the coexistence of both stable and unstable 
Lyapunov exponents in the direction transversal to the 
invariant subspace S Therefore, a closer exami- 

nation of the Lyapunov spectrum could reveal some new 
aspects of the synchronization behavior. We thus con- 
sider the Lyapunov dimension of the system, which is a 
suitable concept to study the Lyapunov spectrum and is 
defined as follows. Let Xj (j = 1, 2, • • • ) denote the j-th 
largest Lyapunov exponent of the system and p be the 



largest integer for which Y^j=i ^ s non negative. Then 
D is given by £| 



D 




|A p+ i 



if there is no such p 
itp = N 



(2) 



In Figure we depict the Lyapunov dimension D versus 
the strength parameter e, for N — 5 and N = 385, and 
for two values of a. We can observe that D assumes a 
maximum D max = N for small values of e. As e enters 
into the synchronization domain, identified in the figure 
by vertical lines, we can observe two distinct behaviors. 
For large N, it monotonically decreases, but for small N 
there is a value e — e^ above which D starts to increase. 
Figure [3] (b) shows the dependence of D with N, for a 
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FIG. 5: Lyapunov dimension D versus e with a = 0.0 and 
a = 0.6. (a) N = 5; (b) N = 385. The vertical lines indicate 
the boundary of the synchronization domain. Inset: detail of 
(b) for 0.8 < e < 1.0 

and e fixed within the synchronization domain. We can 
observe that the Lyapunov dimension tends to saturate 
with increasing N. 

We point now the great similarity between the behav- 
ior of the Lyapunov dimension within the synchroniza- 
tion domain and the behavior observed for the averaged 
synchronization time (t s ), as can be seen by a direct com- 
parison between the shapes of figures [5] and [5] or between 
figures 0(a) and (6). These plots suggests that there is a 
correlation among the Lyapunov dimension and the aver- 
aged synchronization time. In Figure[S]we plot three dis- 
persion diagrams (t s ) x D, each one with two of the three 
parameters a, e, N fixed. All these diagrams give corre- 
lation coefficients p very close to 1, which indicates a very 
strong correlation among these two quantities. In these 
plots the dashed lines correspond to the fitting functions. 
In the first case the fitting is linear and in the two last 
they are exponentials. As the Lyapunov dimension can 
be analytically determined, this result allows also an ana- 
lytical estimation of the synchronization time, which oth- 
erwise requires numerical simulations. The origin of such 
strong correlation between so diversely defined quanti- 
ties is a point that would need a deeper analysis, which 
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is out the scope of this report. Nevertheless, we could 
try to understand this fact on some intuitive grounds by 
observing that it is reasonable to think that the domi- 
nance of negative(positive) Lyapunov exponents in the 
direction transversal to the invariant subspace S would 
tend to minimize (maximize) the synchronization time. 
Accordingly, as it is straightforward from its definition, 
this is precisely the behavior of the Lyapunov dimension. 
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FIG. 6: Dispersion diagrams for (t s ) x D: (a) a = 0.0, £ = 1.0 
and N ranging from N = 5 to N = 2000, p = 0.9998714; (b) 
a = 0.0, N = 501 and e ranging from e = £ c (0, 501) to s= 1.0, 
p = 0. 9849924; (c) JV = 501, £ = 1.0 and a ranging from a = 0.0 
to a=0.2, p = 0. 9868354. 



eraged synchronization time (t s ) and with respect to the 
statistical distribution of the distances d n , which allowed 
us to set N = 385 as a practical lower bound for this sys- 
tem to be considered in the thermodynamic limit in nu- 
merical simulations. We argued that the system behavior 
at the thermodynamic limit can be reasonably well ap- 
proximated, both qualitatively and quantitatively, by its 
behavior at this lower bound. By the way, the behavior 
of systems exhibiting long range couplings, specially con- 
cerning their thermodynamical aspects, is still not well 
understood. We hope that these results could give some 
useful contribution to this subject. 

The Lyapunov dimension for the system within the 
synchronization domain was studied and the results 
showed a very strong correlation among it and the aver- 
aged synchronization time. This result seems interesting 
because it allows an analytical estimation of the synchro- 
nization time, which otherwise requires numerical simu- 
lations. The origin of such a correlation and its related 
consequences are subjects that still need more clarifica- 
tions and it will be postponed to future works. Finally, 
we remark that the statistical distribution of distances in 
Figure |21 can be well fitted by a power law. Further stud- 
ies on this and other related scaling laws on this system 
are now in curse and shall be presented elsewhere. 
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Summarizing, in this report we numerically simulated 
the behavior of a CML consisting of a chain of chaotic 
logistic maps exhibiting power law interactions. We ob- 
served size dependent behaviors with respect to the av- 
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